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ABSTRACT 

Precise near-infrared radial velocimetry enables efficient detection and transit 
verification of low-mass extrasolar planets orbiting M dwarf hosts, which are 
faint for visible-wavelength radial velocity surveys. The TripleSpec Exoplanet 
Discovery Instrument, or TEDI, is the combination of a variable-delay Michelson 
interferometer and a medium-resolution (R=2700) near-infrared spectrograph on 
the Palomar 200" Hale Telescope. We used TEDI to monitor GJ 699, a nearby 
mid-M dwarf, over 11 nights spread across 3 months. Analysis of 106 independent 
observations reveals a root-mean-square precision of less than 37 m s^^ for 5 
minutes of integration time. This performance is within a factor of 2 of our 
expected photon-limited precision. We further decompose the residuals into a 
33 m s~^ white noise component, and a 15 m s~^ systematic noise component, 
which we identify as likely due to contamination by telluric absorption lines. With 
further development this technique holds promise for broad implementation on 
medium-resolution near-infrared spectrographs to search for low-mass exoplanets 
orbiting M dwarfs, and to verify low-mass transit candidates. 

Subject headings: Extrasolar Planets, Data Analysis and Techniques, Astronomical 
Instrumentation, Astronomical Techniques 
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Introduction 



M-dwarf stars are highly compellinja; targ ets for extrasolar planet surveys (iGaidos et al. 



2007 



Lunine et al. 



2009 



Charbonneau 



20091 ) ■ As host stars for transiting exoplanets, M 



dwarfs' low masses, radii and luminosities provide opportunities to broadly characterize 
low-mass exoplanets. Consider that an exoplanet with the same mass, radius and incident 
stellar flux as the Earth orbiting a 0.21 Mq, 0.27 Rq, 3170 K, M5V star generates a 0.69 
m s~^ semi-amplitude radial velocity signal on the host, has a 0.1% transit depth, and 
a 0.01% mid-infrared secondary eclipse depth. These signals are nearly within reach of 
current ground and space-based astronomical instrumentation, enabling measurement of the 
exoplanet's orbit, mass, radius, average density and surface gravity, as well as the planet's 
atmospheric temperature, the role of a greenhouse effect and day/night temperature 
contrasts. Compare this to an Earth-Sun analog system, which produces a 0.09 m s~^ radial 
velocity semi-amplitude, a 0.008% optical transit depth, a 0.0004% mid-IR secondary eclipse 
depth. The lower signals from Earth twins orbiting Sun-like stars limit characterization 
to a measurement of the exoplanet radius with space-based photometry. Current radial 
velocity capabilities and transit timing techniques can put limits on the masses of Earth-like 
exoplanets discovered to transit Sun-like stars. However, accurately measuring their masses 
and surface temperatures will require broad progress in instrumentation. Simply stated. 
Earth-like exoplanets orbiting M dwarfs will be much easier to characterize. 

Unfortunately, the advantages M dwarf hosts provide for characterizing low-mass 
exoplanets come with the challenges of initial detection. M dwarfs are upwards of 5 
magnitudes fainter in My than Sun-like stars, significantly reducing the efficiency of 
visible-wavelength radial velocity planet surveys. Being intrinsically fainter, fewer bright 
M dwarfs occupy a single field on the sky, which reduces the efficiency of single-field 
transit surveys. This is true despite the fact that M dwarfs dominate stellar populations 
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(IWielen et al 



Chabrieiil2003l). Co mbining the My luminosity function for single main-sequence field stars 



19831 ) with the distance modulus, one can calculate the relative number of 
single main-sequence stars in a magnitude-limited wide-field transit survey. Comparing the 
number of single, main sequence stars with My from 2.5 to 8.5, corresponding to F, G and K 
dwarfs, to the number with My from 8.5 to 18.5, corresponding to M-type dwarfs, one finds 
that F, G and K dwarfs outnumber M dwarfs in magnitude-limited wide-field survey by 
more than 1000 to 1. To overcome this limitation, the MEarth t ransit survey is individually 
targe ting bright M dwarfs, rather than staring at a single field (iNutzman fc Charbonneau 



20081 ) at the expense of continuous coverage. 



Despite the challenges, both visible-wavelength radial velocity surveys and transit 
surveys have had success with M dwarfs. As of now, 27 planets have been detected around 
21 M dwarfs, and two are known to transit their host stars: GJ 1214 b and GJ 436 b. All 



hosts have a spectral type earlier than M5. GJ 1214 b was initial 



transit survey, then confirmed by radial velocity measurements (ICharbonneau et al. 



ly detected by the MEart h 



and GJ 436 b was initially discovered by a radial velocity s urvey fiButler et al 



20091), 



with the transit detection occurring later ( iGillon et al. 



2OO4J ) 



20071 ). At this stage it is unclear 



whether the majority of future M dwarf transiting planets will be detected first by transit 
surveys then followed up with radial velocity measurements, or vice-versa, given that both 
individually target M dwarfs. However, the Kepler Mission is a sensitive wide-field transit 
survey observing approximately 2500 M dwarfs brighter than V = 14. If low-mass planets 
are common around M dwarfs, Kepler will potentially detect dozens of transiting candidates 
around M dwarfs which are bright enough for follow up radial velocity measurements. 
Recent statistical analysis of Kepler planet candidates indicates that short period planets 



are, i n fact, more common around low- mass dwarfs than high mass dwarfs (IHoward et al. 



201ll ). although this analysis did not include stars later than MO. Measuring precise 



radial velocities of dozens of faint M dwarf hosts will prove a difficult challenge. The 



- 6 - 



effort to detect and confirm more transiting terrestrial exoplanets can be accelerated by 
improving radial velocity precision at near-infrared wavelengtlis. M dwarf spectral energy 
distributions peak in the near infrared, and the many absorption lines from molecular 
transitions in ne ar-infrared M dw arf spectra provide rich structure for measuring radial 



velocity signals (jJones et al 



20091 ). Precise near-infrared radial velocities will increase the 
efficiency of exoplanet detection and transit verification around early M dwarfs, and will 
enable detection and transit verification of exoplanets orbiting M dwarfs later than M5, 
which have yet to be discovered. 

Precise radial velocity measurements of M dwarfs must contend with activity-related 
radial velocity "jitter," int roduced by temperatur e inhomogeneities such as spots or flares on 



the rotating photosphere (ISaar fc Donahue 



19971 ). G- and K-type stars show a co rrelation 



between radial velocity jitter and Ca II emission, an indicator of surface activity (jWright 



20051 ). Ha emission is a typical indicato r of surface activity on M dwarfs, and has been 



shown to correlate with Ca II emission (Walkowicz &: Hawlev 



r 

dwarfs in the Sloan Digital Sky Survey, iWest et al. 



20091 ). Using spectra of M 



( I2OO4I ) found the fraction of M dwarfs 



with Ha emission rises st eeply with 



over 70% for M8 dwarfs flWest et al. 



ater s pectral type, from less than 5% for MO dwarfs, to 



20041 ). This could present a challenge for precise radial 



velocimetry of mid-to-late M dwarfs. 

However, less than 10 m s~l of visible- wavelength radi al velocity performance ha s 



been achieved on early M dwarfs which show Ha emission flEndl et al. 



2003 



Reiners 



20091 ) . And, radial velocity jitter is expected to be significantly reduced at near-infrared 



waveleng ths, because tempe rature inhomogenieties on the stellar surface have lower flux 



contrast (IBarnes et al. 



20101 ). For this reason, precise near- infrared radial velocities also 



provide an import ant tool for verifying planetary candidates identified around active stars 



(jPrato et al. 



20081 ). For example, visible- wavelength radial velocity measurements of the 
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young, active star TW Hydra indicated the presen ce of a planet (Setiawan et al.ll2008l ). 



but infrared radial velocity measurements did not (iHuelamo et al. 



20081 ). leading to the 



conclusion that the visible-wavelength Doppler signal was likely due to rotating spots. 

Measuring precise near-infrared radial velocities is challenging because of the lack of 
available high resolution spectrographs, the difficulty of calibration, and telluric interference. 
Several groups have had success using the Earth's atmosphe ric lines to calibrat e high 



Blake et al. 



teoioh achieved 



resolution spectrographs, similar to the iodine cell technique: 
50 m s^^ of precision on an M dwarf over several years with th e NIRSPEC spectrograph 



(R=25000) on the Keck II telescope, and 



Figueira et al. 



(I2OIOI ) achieved 6ms ^ on an 



Sun-like star over one week with the CRIRES spectrograph (R= 100000) on the European 
Southern Observatory's Very Large Telescope. So far, the best precis i on ha s been achieved 



Bean et al. 



( I2OIOI ) achieved 5 m 



by calibrating CRIRES with an ammonia gas cell, where 
s~^ of long term precision on an M dwarf. 

High resolution near-infrared spectroscopy is becoming an important technique for M 
dwarf exoplanet science, but in order to significantly impact the field it must be competitive 
with visible-wavelength spectroscopy on M dwarfs. Consider GJ 1214, a V=15.1 M4.5 



dwarf on which HARPS achieved 5ms ^ o: 



minute integration time (ICharbonneau et al. 



la radial velocity uncertainty with a 40 



2OO9I). Scaling the CRIRES resuh to GJ 1214, 



the ammonia-cell technique would achieve 15 m/s of Icr radial velocity uncertainty with 
the same exposure time as HARPS. This is because current high-resolution near-infrared 
spectrographs suffer from limited simultaneous bandwidth, with NIRSPEC covering 400 
nm at one time and CRIRES covering 80 nm. To overcome this, several groups have 
proposed or are developing high-resolution near-infrared echelle spectrographs with large 



simultaneous bandwidth, such as the Pre cision Radial Velocity Spectrograph f|Jones et al. 



20081 ) . the Habitable-Zone Planet Finder (IMahadevan et al 



201G 



Ramsey et al. 



20081), and 
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CARMENES flQuirrenbach et al.lboioh . 



High precision velocimetry and large simultaneous bandwidth can also be obtained 
by introducing a Michelson interferometer into the optical path between a telescope and 



an existing me dium-reso . 



interferometrv (lErskine 



van 



2007 



vken et al 



Behretal 



ution spectrograph, a method known as ext ernally disperse d 



20031 ). dispersed fixed delay interferometry (IGe et al. 



2010 ) and also dispersed Fourier transform spectroscopy (IHajian et al 



2006 



20091 ). The interferometer multiplies the stellar spectrum by a sinusoidal 



transmission comb before being dispersed by the spectrograph, which creates a moire fringe 
pattern highly sensitive to Doppler shifts. Medium-resolution near-i nfrared spectrograp hs 



are becoming 



Simcoe et al. 



wide ly available on 4-to-6 meter-class telescopes (e.g. 



Wilson et al. 



2004 



20101 ) because of their smaller size and complexity compared to high-resolution 



spectrographs. The availability of medium-resolution near-infrared spectrographs makes 
externally dispersed interferometry a promising technique for broad implementation to 
search for these most interesting low-mass exoplanets. 

In this paper we present results from the TripleSpec Exoplanet Discovery Instrument, or 
TEDI, designed to measure precise near-infrared radial velocities of nearby M dwarfs. TEDI 
is the combination of a Michelson interferometer and TripleSpec, a facility near-infrared 
spectrograph on the Palomar 200" Hale Telescope that simultaneously covers 1.00 to 2.46 
fim at a resolution of 2700. 



2. Theory 



The theory behind externally dispersed interferomet rv or disperse c 



spectroscopy is c 


escri ; 


van Evken et al. 


2010) 



Ded in several previous papers (e.g. lErskine 



2003 



Fourier transform 



Hajian et al 



2007 



20101 ). These treatments are often specific to the instruments involved. 
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and do not account for some of the effects that we encounter with TEDI. Here we provide a 
theoretical accounting of the TEDI data product and an explanation of how that product 
is converted into a measured change in radial velocity. 



2.1. The TEDI Data Product 

Given wavenumber u = 1/A, an intrinsic stellar spectrum S^-, spectrograph line-spread 
function Ry, and an interferometer optical path difference of r, an individual EDI spectrum, 
ly^T) can be described as: 



= [Su{l + cos(27rrz/))] * (1) 

where * indicates a convolution and (1 + cos(27rrz/)) represents the effect of the sinusoidal 
transmission comb introduced by the interferometer before convolution with the spectrograph 
line-spread function. The spectrograph line-spread function represents the broadening of 
spectral lines introduced by the spectrograph. In order to get moire fringes, the delay 
T must be slightly varied by At to modulate the resultant spectra. This can be done 
by either positioning one of the interferometer mirrors such that different delays appear 



perpe ndicular to the dispersion direction of the spectrograph (e.g. 



Erskin 



2003 



Zhao et al. 



1 

2009.), or by actively moving one of the interferometer mirrors with a piezo actuator and 



taking an individual spectrum at each position. With TEDI, we use the latter approach, 
moving the mirror to keep the interferometer modulation within a pixel rather than across 
several pixels. This dramatically reduces the effects of pixel-to-pixel calibration errors, such 
as effects from poor flat-fielding and/or background subtraction. 

If we rewrite r as a bulk delay plus a small phase shift (tq + At) and assume that 1 / At 
is large compared to a resolution element of the spectrograph, then we can remove At from 



the convolution integral. In TEDI, l/Ar is 3 orders of magnitude larger than a resolution 
element. Applying this and rearranging Equation 1 using trigonometric identities, we can 
express an individual TEDI spectrum as a function of the bulk delay and phase shift: 



L,ToAr = + By cos(27rAri/) - sin(27rAri/) (2) 

where: 

A,^S,* R, (3) 

= [5,cos(27rToi/)] (4) 

a = [>5,sin(27rToi/)] (5) 

Ai, is the spectrum of the star at the native-resolution of the spectrograph, referred to as the 
"conventional" spectrum. and describe the moire fringes that contain the signal from 
the high-resolution stellar features. By slightly changing the delay of the interferometer by 
steps of Ar, and taking spectra at each step, we can fit A^, B^, and Cj, to the modulation 
at each pixel according to Equation 2. We can then construct the complex visibility: 

B, - iC, = [S.e-'^^^'''] * K (6) 

The complex visibility contains a real and imaginary component which can also be described 
by a phase, = arctan (Bj^/Cj,), and a visibility, K = ^/B^ + C^- Under a small Doppler 
shift of the stellar spectrum Siy+Au, where and assuming that Au 

and 1/tq are small compared to a resolution element of the spectrograph, Au can be 
transferred to the exponent by treating the convolution as a Fourier transform integration 
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and applying the Fourier shift theorem. In this case, the Doppler-shifted epoch complex 
visibility, Bl — iC^, is related to an unshifted template complex visibility, -B° — iC^, by: 

Bl - iCl = [5° - .COje-^^TrroA. (7) 

Thus, a small Doppler shift causes a change of the phase of the complex visibility. Since all 
wavelengths expect the same radial velocity shift, the change in phase versus wavelength is 
expected to follow a simple curve: 

A(/) = 27rToAi/ = 27rTo(ARV)i//c (8) 

where ARV is the change in the radial velocity of the star. However, the motion of the 
telescope relative to the barycenter of the Solar System introduces large Doppler shifts on 
the order of 10 km s~^ on most stars, which is a significant fraction of a TripleSpec pixel. 
In this case, the convolution cannot be treated as Fourier transform integration and the 
Fourier shift theorem is not apphcable. To accurately account for a large change in the 
radial velocity of a star, we define a complex line-spread function as Ri, = e~^^'^'^°^^^R^. We 
can now rewrite Equation 6 as: 

B, - iC, = e^2.ro.[^^ ^ 

Shifting a function before a convolution is equivalent to shifting after a convolution; 
therefore, a Doppler shift in Si, will cause the quantity [Si, * R^,] to shift by the same 
amount. The Doppler shifted epoch complex visibility, Bl + iCl relates to an unshifted 
template complex visibility, B^ + iC^, by: 
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Bl - iCl = e^2-o-[(50 _ ^dy-^'^'^-^X-ru^t^^ (10) 

The change in radial velocity, ARV, between the template and epoch can be measured 
by multiplying the template S° — by e"'^'^'^"'', interpolating the product onto a 
Doppler-shifted v + Ai/ grid, multiplying that by e*^'^'^"'', and then comparing the resulting 
complex visibility to an epoch complex visibihty {B\, — iCl). The visibihties, K, of the 
template and epoch measurements depend on the quality of interference or contrast in the 
interferometer in addition to the radial velocity shift, but the phases (p^ only depend on the 
radial velocity shift. For this reason, we fit a shifted template to an epoch measurement with 
ARV as a free parameter, using the error-weighted difference in phases as the goodness-of-fit 
statistic. 

In configuring TEDI for a particular target, the value of the bulk delay Tq is chosen to 
produce the largest radial velocity signal in Bi, and C^,. This depends on several factors 
with the rotational broadening of the target's absorption lines playing the largest role. The 
rotational broadening relates to the rotation of the star, Vrot, projected along the angle of 
the rotation axis with respect to the line of sight, i, such that Az/jot = (Vrot sin i)^. The 
radial velocity signal is maximized when the periodicity of the interferometer comb, I/tq, 
matches the width of the stellar lines, Ai/rot, causing the largest modulation of spectra, 
Iu,to,At-^ with At. Rotational broadening increases linearly with wavenumber, but the period 
of the comb does not. This means that each wavelength has a different optimal bulk delay; 
optimizing the bulk delay across an entire spectrum is typically done by choosing the bulk 
delay corresponding to the spectral region with the highest line density and highest flux, 
both of which depend on the spectral type. For these reasons, the bulk delays available 
in TEDI were chosen based on performance simulations using rotationally broadened high 
resolution models of late-type stars provided by Travis Barman using the PHOENIX model 
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[e.g. 



Fuhrmeister et al 



20051 ) ■ We used a bulk delay of 4.6 cm for the results in this paper, 
corresponding to a projected rotational velocity of 1 to 5 km s^^ for a mid-M star. 



2.2. Determining the Interferometer Delay 

In our account of the TEDI data product and ARV measurement, we have relied 
on accurate knowledge of the bulk delay and phase steps, which together make up the 
interferometer delay. However, errors in these quantities will result in errors in the measured 
change in radial velocity. An error in the estimation of the bulk delay will appear as a 
proportional error in ARV. For example, a 1% error in Tq will correspond to a 1% error 
in ARV. These errors can be calibrated by observing radial velocity standard stars, and 
changing the bulk delay to minimize the residuals between the ARV measurements and 
those expected from the motion of the telescope relative to the Solar System barycenter. 

An error in the estimation of the phase steps is much more severe. If all of the phase 
steps used to construct the complex visibility are incorrectly estimated by an offset, 5(Ar), 
this will correspond to an offset in the radial velocity change, 5(ARV), of 



5(ARV) = (11) 

To 

With a bulk delay of 4.6 cm, a 1 nm offset in the phase steps corresponds to a 7 m s~^ error 
in the measured change in radial velocity. 

Ensuring that the interferometer has a specific delay with 1 nm of accuracy for every 
epoch observation is extremely challenging. Instead, with TEDI we calibrate the bulk 
delay and phase steps with emission lines from a TliAr hollow cathode lamp, and use the 
modulation of the emission lines to correct for delay differences after the data has been 
collected. 
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2.3. Telluric Calibration 

Spectroscopic observations at near-infrared wavelengths must contend with telluric 
absorption lines introduced by the Earth's atmosphere. Telluric lines are caused by many 
species of molecules in the Earth's atmosphere, and the largest contributors at near-infrared 
wavelengths are H2O, CO2 and CH4. Telluric lines are numerous across near-infrared 
wavelengths and will contaminate the stellar complex visibility. In conventional spectroscopy 
the standard method for removing telluric lines is to observe a featureless star with a 



known spectrum, such as a rapidly rotating B star or an AO, Vega- 



Vacca et al 



ike star, at a similar 



2003f ). Unfortunately, 



airmass, and divide that into the target spectrum (e.g. 
the complex visibility is not amenable to a simple division, necessitating a more complex 
and less robust calibration technique. Simulations of telluric contamination and removal 
indicate the best method to empirically calibrate telluric lines in TEDl observations is 
to subtract the complex visibility of a calibrator from that of the target, normalized to 
the conventional spectrum Ai, of the target. However, incomplete or inaccurate telluric 
calibration will introduce an error into the measured radial velocity. Simulations indicate 
that telluric calibration currently limits TEDl performance, which we address in Section [6l 



TEDI 



A thorough 



(lEdelstein et al. 



description of the TEDI design and hardware is available elsewhere 



20101), so here we briefly summarize only those portions of the instrument's 



design, components, and history relevant to the cu rrent performance. The TripleSp ec 



2004 



Herter et al. 



20081). 



spectrograph was commissioned in October of 2007 (1 Wilson et al. 
TripleSpec is a cross-dispersed long-slit near-infrared echelle spectrograph, dispersing a 1 x 
30 arcsecond slit from 1.0 to 2.5 across 5 orders at resolution 2700 onto 2 quadrants of 
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a Rockwell Scientifico Hawaii-2 HgCdTe detector. 

In December of 2007, we attached an interferometer to TripleSpec to form TEDI. The 
original design used a mirror to redirect the telescope beam toward the interferometer, 
with a pellicle beam splitter used to simultaneously inject TliAr emission light. The design 
su ffered from non- c ommo n path errors between the starlight and ThAr light, described 



m 



Muirhead et al. 



( 120101 ). In December of 2009, we removed the interferometer from 
TripleSpec and upgraded the design to eliminate the non-common path behavior. The 
overlapping starlight and ThAr light is now focused onto a fiber before the interferometer, 
which ensures common path and common delay between the two. The new design includes 
two fibers; nodding the target between the fibers nods the target on the TripleSpec detector, 
allowing for efficient background subtraction. The upgraded interferometer was attached to 
TripleSpec in June 2010. 



3.1. The TEDI Beam Path 

Figure [1] depicts the beam path through the interferometer before entering TripleSpec. 
First we describe the injection of starlight and ThAr light into the science fibers: (1) A 
mirror mounted on a removable swing arm defiects the f/16 telescope beam that would 
otherwise go to TripleSpec. (2) A mirror mounted on a 3-axis, tip/tilt/piston piezo actuator 
directs the beam towards the two science fibers. (3) A visible/near- infrared dichroic refiects 
the visible light to the CCD camera (4) and transmits the near-infrared light to the science 
fibers. Once a star has been placed on a science fiber, the CCD camera (4) is used to 
monitor the location of the star's visible image, and tip/tilt commands are sent to the piezo 
actuator in (2), ensuring that the star does not drift off of the science fiber. To introduce 



Now Teledyne Technologies 
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the ThAr calibration light a small telescope (5) images a 1 mm fiber core carrying the ThAr 
emission fight to the science fibers at f/16, matched to the telescope beam. The magnified 4 
mm fiber core illuminates both science fibers with ThAr light. (6) A fiat window mixes the 
ThAr fight into the science beam by refiection. (7) The science fibers accept the starfight 
and ThAr fight via glued-on microlenses, which convert the f/16 beam to an f/4 beam. 

Next we describe the interferometer: (8) Starfight and ThAr fight exit the fibers 
via an identical set of microlenses at f/16, with an additional single- mode fiber carrying 
He-Ne laser light mounted nearby. (9) An ofi^-axis parabola coUimates the beams. (10) 
The beam splitter directs half of the beam to a fixed mirror, and the other half through 
one of several etalons (11) to a mirror mounted on a second 3-axis piezo actuator and an 
adjustable linear stage (12) used in combination with the etalon to introduce the delay. 
(13) One interferometer output is directed to another visible/near- infrared dichroic (14), 
which sends the near-infrared light to a chopper/photodiode system to actively monitor 
throughput (15). The second dichroic sends visible light to a second CCD camera (16). 
The second CCD camera images the interference pattern produced by the He-Ne laser light. 
The pattern is altered by adjusting the tip/tilt /piston piezo actuator in (12). When the 
interference pattern contains no fringes, the interferometer is well aligned and produces the 
highest stellar moire visibifities. The other interferometer output is directed back toward 
the TripleSpec window via a fold mirror (17), with an adjustable mirror located at the 
pupil (18) for positioning the image of the fiber on the TripleSpec slit. A mirror underneath 
the removable swing arm (19) redirects the beam down into TripleSpec at f/16. The entire 
interferometer assembly sits above TripleSpec at the Cassegrain focus and fits inside the 
39" diameter hole in the 200" primary mirror. 
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3.2. Throughput 

TEDI adds many optical surfaces into the beam path prior to entering TripleSpec, 26 
surfaces for one interferometer arm and 30 for the other, and uses 50% of the starhght for 
monitoring throughput. This is necessitated by the geometry of adding the interferometer to 
TripleSpec at its existing Cassegrain mount, but results in significant throughput loss. With 
1 arcsecond seeing, the total throughput of the system-including the atmosphere, telescope, 
interferometer and TripleSpec-peaks at 1.5% at 1.68 /im. Without the interferometer, the 
throughput peaks at 30%. Focal ratio degradation in the optical fibers also contributes to 
the lower throughput. TripleSpec contains a Lyot stop which, when combined with the slit, 
matches the seeing-limited etendue of the telescope. Focal ratio degradation in the fibers 
will effectively increase the etendue of the beam, and introduce irrevocable losses at either 
the Lyot stop or the slit, depending on the focus. 



3.3. Fiber Scrambling 



The output illumination of a multimode fiber is dependent on various factors 
including the input illumination geometry. Fiber scrambling methods can be applied to 
reduce this dependence so that the output fiber illumination is less affected by changes 
in seeing, telescope pointing or tel e scope focus . Met hods include an optical double 



scrambler (iHunter &: 



( iBaudrand fc Walkei 



Ramsey 



2001 



19921: 



Lovis et a. 



Ramsey et al 



20061 ) or a mechanical fiber agitator 



20081 ). though the methods are not equivalent 



and produce different fiber output illumination. TEDI does not use supplemental fiber 
scrambling methods. 

Partial fiber scrambling is achieved by the microlenses attached to the science fibers, 
which image the telescope pupil onto the fiber core and act as an image pupil exchanger. 
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similar to a "double-scrambler" ( iHunter &: Ramsey 



1992I ). However, without identical input 



illumination, the ThAr and stellar cavity illumination will behave slightly differently. The 
ThAr injection is mechanically fixed while the stellar injection will change depending on 
seeing, guiding and focus fluctuations. This could introduce errors into the calibration of tq 
and At, if the slightly different illumination through the interferometer results in different 
optical paths. We measured this effect experimentally by simulating fiber illumination 
fluctuations and believe it is not a significant source of radial velocity fluctuations. We 
discuss the experiment and results in Section 16.11 



4. Observational Procedure and Data Analysis 

Figure [2] shows a sample TEDI exposure with both mixed star-ThAr light and ThAr 
light alone. A single TEDI measurement of the complex visibility in one fiber consists of 20 
spectra: 10 of mixed star-ThAr light, and 10 of ThAr alone. After 10 exposures with the 
star focused on one fiber, the star is nodded to a second fiber. Since both fibers receive 
constant ThAr light, the nodding procedure efficiently provides 10 mixed star-ThAr and 10 
ThAr alone spectra on both fibers. Between each of the 10 exposures we change the delay 
by 0.25 fim using a closed- loop Piezo actuator (element #12 in Figured]). In order to fit 
the coefficients A^,, B^,, and C^, in Equation 2, At must Nyquist sample Iu,ToAr f*^^ ^■ 
Ten spectra with steps of 0.25 fim ensures this condition for the full TEDI bandwidth., and 
we refer to this set of 10 spectra as a phase set. 

Converting 20 exposures, and 40 spectra, into complex visibilities for both fibers 
involves several steps: spectral extraction, deconvolution of the ThAr lines from the mixed 
star-ThAr spectra, fitting A^,, B^, and Ci, to both the ThAr and stellar spectra, accounting 
for changes in the start delay between template and epoch by referencing A,^, By and 
to the same start delay, and subtracting the complex visibility of a telluric calibration star. 
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With the complex visibihties measured for a template and an epoch, we measure the radial 
velocity change between the two by fitting a shift as in Equation 9. We treat each fiber 
independently, with independent template and epoch measurements. 



4.1. Spectral Extraction 

Each TripleSpec exposure contains several electronic anomalies which were removed 
during data processing by modeling the effects and subtracting them out. An IDL routine 



for removing the electronic anomalies in TripleSpec exposures is availab 



website 



and this routine has been used in pre vious TripleSpec results ( jMiller et al. 



e on the author's 



2011 



Rojas-Ayala et al. 



2010 



Martinache et al 



20091 ). The details of the electronic anomalies are 



beyond the scope of this paper, but they involve crosstalk of signals between the detector 
quadrants, and capacitive coupling of signals between the channels within each quadrant. 
For a discussion of crosstalk and capacitive coupling in nondestruc tive ar r ays, i ncluding 
Hawaii-2 HgCdTe detectors, we refer the reader to Chapter 6.6.4 of iRiekd ( 12002| ). 



The TripleSpec detector is subject to a variable electronic background, which changes 
slowly across the detector and slowly with time. Th e electronic background is mostly 



removed by subtracting multiple correlated readouts (IGarnett fc Forrest 



19931 ). However, 



correlated sampling does not fully remove the electronic background in the detector, which 
persists and changes with each exposure. To remove the fluctuating electronic background, 
we subtract the median of a 20 x 20 pixel box centered around the extraction point. Each 



exposure was flltered for cosmic rays and hot pixe 



available from the IDL Astronomy User's Library (iLandsman 



s using the IDL routine sigma_fllter. 



19931 ). which locates pixels 



signiflcantly brighter than their neighbors. We flag these pixels as bad, and they are not 



http : / /astrosiin2 . astro . Cornell . edu/~muirhead/tspec_clean . pro 
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used in any further analysis. 

Individual exposures are flat-fielded and dark subtracted using dome flats and dome 
darks taken before an observing night. It is difficult to attain a true estimate of the read 
noise because TriplcSpec does not have the capability to take a "cold dark," wherein a 
cryogenic mask fully occults the detector. Instead, we use the dome darks to construct an 
image of the read noise of the detector, as well as flag any hot or variable pixels. Dome 
darks contain significant background emission in K band, and subtracting them from target 
exposures leaves a negative K band background. We remove the negative background after 
the spectra are extracted, however the photon noise from the thermal emission in the dome 
dark inflates the estimate of the read noise. 

We extract the spectra by (non-optimally) summing pixels across the slit image on 
the TripleSpec detector. The TripleSpec slit is tilted on the detector by different amounts 
at different wavelengths, and summing pixels in this fashion slightly degrades spectral 
resolution. The tilt is nowhere more than 10 degrees; as such the resolution loss is less 
than 3%. We chose to sum pixels along the slit to avoid aliasing effects introduced when 
interpolating the slit onto a rectilinear grid, or interpolating a profile onto the non-rectilinear 
detector. Wc replace bad pixels with the value of a profile fit to the slit image, where the 
profile is constructed from nearby slit images. For slit images where more than 40% of 
the pixels are flagged, the entire wavelength channel is flagged as bad. For each extracted 
wavelength channel, we subtract the median of a 20 x 20 box centered on those pixels to 
eliminate a fluctuating electronic background in the TripleSpec detector. 

A two-dimensional polynomial wavelength solution for the TripleSpec detector was 
determined using gas discharge lamps and a slit mask during commissioning of the 
spectrograph. TripleSpec has no internal moving parts, which minimizes flexure and 
changes in the wavelength solution. However, large changes in the gravity vector of 
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TripleSpec result in illumination shifts of approximately 1 pixel. After extraction, Gaussian 
profiles are fit to bright, isolated ThAr lines, and the profile centroids are used to offset to 
the original wavelength solution to account for flexure. Figure E] plots an example mixed 
star-ThAr TEDI spectra, in units of signal-to-noise. 

In order to measure radial velocity changes, it is more important that the wavelength 
solutions of the spectra are the same for every measurement, rather than accurate for every 
measurement. This is ensured by cross-correlating the A^^s of the ThAr-alone spectra of the 
epoch measurements onto the A^, of the ThAr-alone spectra of the template measurement. 
The shift applied to the epoch ThAr spectra is then applied to the corresponding mixed 
Star-ThAr epoch spectra. For sub-pixel shifts, the spectra are shifted using spline 
interpolation. 



4.2. Star/ThAr Separation 

To measure the set of phase steps introduced into a mixed star/ThAr phase set, 
the ThAr lines must be separated from the mixed star/ThAr spectra. This is done by 
interpolating the stellar spectrum underneath the ThAr lines and then subtracting that 
interpolated spectrum. The locations of the ThAr lines are determined using the ThAr-alone 
phase set taken immediately before or after the mixed spectra on the same pixels using 
the nodding scheme. From here on, we make the distinction between ThAr-alone and 
ThAr-separated phase sets, the former being the phase set of ThAr alone immediately 
preceding or following the mixed star-ThAr phase set, and the latter being separated from 
the mixed star-ThAr spectra. Figure H] plots images of the 10 mixed and separated stellar 
and ThAr spectra, and shows the resultant moire fringes, to which Ai,, and are fit. 



- 22 - 



4.3. Fitting A^, and 

In order to fit the coefficients A^j, By and Cp to the stellar spectra, we first fit them to 
the ThAr-alone and ThAr-separated phase sets, which determines the sizes of the phase 
steps. At. The model for the data is Equation 2, with A,^, 5,^, and 9 of the 10 phase 
steps, Ar, as the free parameters. The first phase step is fixed at 0. The fitting procedure 
is a joint non-linear Levenberg-Marquardt (LM) fit for the phase steps and a linear fit for 
Ay, By and C^, at each wavelength channel. For each LM iteration of the 9 phase steps, a 
linear fit of Aj,, B^, and Cp is performed for each wavelength channel of the phase set. The 
best-fit delays, and best fit Ay, By and Cy, are found when the error- weighted residuals 
between the model and data are minimized. The fitting procedure adjusts the phase steps 
to ensure that the ThAr lines fiuctuate sinusoidally. This method was chosen over a fully 
non-linear fit to increase computational efficiency, since there are over 10000 wavelength 
channels in a given phase set. A gen eral description of the joint LM-linear fitting method is 



available in iWright fc Howard! (j2009|) . 



When fitting Ay, By and Cy to the stellar spectra, the phase steps are fixed since 
they are determined by the separated ThAr lines. Therefore, only the linear portion of the 
fitting routine is used. Before fitting the stellar coefficients, the ThAr-alone coefficients are 
used to create a model of the background at the phase steps of the mixed star-ThAr phase 
set, and that is subtracted from the mixed star-ThAr spectra. This removes the thermal 
background and OH airglow lines present in the spectra, and is equivalent to subtracting a 
nodded pair of exposures while still accounting for the effect of the interferometer. 

The stellar spectra ffuctuate due to changes in seeing and telescope guiding errors, 
which introduce ffux variations that are a function of wavelength. To correct for this, the 
average shape of each of the 10 spectra is found by convolving each spectrum with a running 
boxcar mean filter. The shapes are used to normalize each of the ten spectra and remove 
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these fluctuations. Various widths for the fllter were tested to find that which produced the 
lowest radial velocity residuals. If the width of the filter is too large, it will not correct for 
fast fluctuations versus wavelength, and if it is too narrow it will remove the modulation 
introduced by the interferometer, which is the signal itself. We found the width with the 
lowest residuals to be 11 wavelength channels, or 4 resolution elements. After the stellar 
spectra are normalized, A^,, and are linearly fit to each wavelength channel. 

4.4. Referencing Delays 

The LM-linear fitting routine determines the sizes of the phase steps between each 
spectrum, by ensuring sinusoidal variation of the ThAr lines, but does not determine the 
difference between the start delay of one set and the start delay of another set. To reference 
a complex visibility to a different start delay, we multiply it by an exponential phase change: 

[B, + iC,]^[B, + iC^y'''' (12) 

where 5t is the change in delay. To correct for the difference in start delays between 
two complex visibility measurements, we fit for the 5t which reduces the phase difference 
between the separated ThAr complex visibilities of a template and epoch measurement. 
That 5t is then applied to the epoch stellar complex visibility, such that the template and 
epoch are referenced to the same start delay. 

4.5. Telluric Calibration 



Telluric calibration is performed by measuring the complex visibility of a rapidly 
rotating A-type or earlier "standard" star near the target, where only the telluric fines will 
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contribute to the complex visibility, and subtracting the telluric complex visibility from 
that of the target. The standard complex visibility must be normalized to the target to 
account for the differences in flux, interferometer visibility and airmass, which will change 
the visibility as the telluric lines change in depth. 

To normalize the standard star's complex visibility, it is first referenced to the same 
delay as the target, then divided by of the standard and multiplied by A^, of the target, 
which accounts for the change in flux. Then, a boxcar filter versus wavelength fits a running 
multiplicative offset to the standard complex visibility to best match the target complex 
visibility for a given wavelength bin. A boxcar width of 51 pixels was found to give the 
lowest radial velocity residuals. 



4.6. Measuring Radial Velocity Changes 

As mentioned previously, a radial velocity change is measured between a template and 
an epoch measurement. The template measurement is identical to an epoch measurement, 
and is simply chosen as a radial velocity zero-point. With template and epoch A^, B^, and 
C*!/, all referenced to the same delay start point, and the telluric effects at least partially 
removed, the radial velocity difference ARV can be measured using the procedure outlined 
in Section [2?T1 This requires knowledge of the bulk delay tq. In order to measure the bulk 
delay, the phase of the complex visibility of the ThAr spectra must be modeled, which 
requires accurate knowledge of the wavelengths of the ThAr lines. Using a high-resolution 
spectrum of a ThA r emission lamp obtained for calibrating the CRIRES spectrograph 



(iKerber et al. 



20081 ) ■ we directly model the phase of By and Cy using several well-separated 
ThAr lines in J band, with the bulk delay as the only free parameter. The delay of the 
ThAr alone set for the template, which all subsequent sets are referenced to, is chosen for 



calculating the bulk delay. As stated in Section 12. 2[ any residual errors in the bulk delay 
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can be corrected using radial velocity standard stars. 

With To known, the template complex visibility is resampled onto a finer wavelength 
grid and multiplied by e"*^'^'^"'", then shifted by Az/ = ^^i^ using spline interpolation. The 
resulting complex visibility is multiplied by e*^'^'^°'", returning it to a slowly varying function 
of wavelength. This is then resampled onto the original wavelength grid, and the phase (p^, 
is compared to an epoch measurement. The ARV which minimizes the of the phase 
difference between template and epoch, weighted by the formal uncertainty in the phase 

is recorded as the measured change in radial velocity. 

Instead of using the full TripleSpec bandwidth of 1.00 to 2.46 fim for measuring ARV, 
we currently limit our measurements to 1.48 to 2.15 fim, ignoring regions which have strong 
telluric absorption. This bandpass has the highest instrument throughput, strong stellar 
absorption lines, and was found empirically to deliver lower residuals than when including 
the full bandwidth. Figure [5] plots the difference in phase converted to radial velocity 
between a ARV-shifted template complex visibility and a typical epoch complex visibility 
of GJ 699 for each wavelength channel in this bandwidth. 



4.7. Formal Uncertainties 

Formal uncertainties for each ARV measurement are found by carrying the per pixel 
errors in each exposure through the reduction process. The error in each wavelength channel 
of each spectrum is calculated by summing in quadrature the read noise and photon noise 
in each contributing pixel. Formal uncertainties in A^, By and are calculated during the 
linear portion of the combined linear-LM fitting routine. The uncertainties in By and Cy 
are passed to the ARV fitting routine and included to weight the phase difference between 
the shifted template and epoch. The formal uncertainty in ARV, o"arv, is calculated by 
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converting the formal uncertainty in the phase of the complex visibility to an uncertainty 
in the radial velocity, and summing the radial velocity uncertainties in each wavelength 
channel in inverse quadrature: 



C^ARV= J^^^ (14) 
V ^ -inv. 

Figure [6] plots a histogram of the residuals and the cumulative distribution function 
of the radial velocities in Figure |5l normalized to their la uncertainties. The cumulative 
distribution function indicates significant outliers. To reduce the effect of wavelength 
channel outliers, we remove those wavelength channels with the highest 3% of residuals 
and repeat the ARV fit. Currently, errors in the estimates of At and errors in the delay 
referencing of complex visibilities are not included in the calculation of ctarv- Errors 
introduced by these effects are presumed to be small, given the high signal-to-noise of the 
TliAr calibration lines compared to that of starlight. 



5. Observations 



To test the precision of radial velocities measured by TEDI, we conducted a 
commissioning campaign o n GJ 699, also known as Barnard's Star. GJ 699 is a bright 



f J=5.24, H=4.83, K= 4.52 



Skrutskie et al 



Browning et al 



20061 ). slowly rotating (Vrotsini < 2.5 km s ^; 



20101) M4 dwarf. Visible-wavelength radi al velocity measu rements limit 



potential radial velocity variations to less than 7.2 m s ^ (lEndl et al 



20031 ). making it an 



ideal candidate for testing the radial velocity precision of TEDI. 
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Each of the 20 exposures which make up two ARV measurements of GJ 699, one for 
each fiber, was 30 seconds long, achieving a median signal-to-noise per wavelength channel 
of roughly 250 in H band for each spectrum. The 20 exposures took approximately 15 
minutes of total observing time including exposure time, read time, time used to nod the 
star between fibers, and time between exposures used to change the delay and communicate 
between the TEDI control computer and TripleSpec control computer. We used the highest 
bulk delay available in TEDI of 4.6 cm, optimized for low projected rotational velocity. 

We took 53 such observations of GJ 699, spread over 11 nights in June, July, August 
and September of 2010. That is, 53 ARV measurements with Fiber A, and 53 with Fiber 
B over 4 observing runs. For telluric calibration, we observed 7 Ophiuchi, a V=3.75 AOV 
star 3.17 degrees a way from GJ 699, with a projected rotational velocity of 210 km s~^ 



( iRoyer et al 



20071 1. At such a high projected rotational velocity, any spectral features of the 
star are significantly broader than the interferometer comb spacing and do not contribute 
to the measured moire pattern and complex visibility. We observed 7 Ophiuchi roughly 
once for every 4 measurements of GJ 699, or once every hour. 

Initially, the ARV residuals correlated proportionally with the expected ARV from the 
motion of the telescope relative to the barycenter of the Solar System. This is a consequence 
of applying an incorrect value for the bulk delay tq in the data reduction procedure, as 
described in Section 12.21 To compensate we altered the tq used in the reduction to minimize 
the ARV residuals, and found a best value of 4.654 cm. For future measurements, ARV 
measurements of known radial velocity standard stars will be used to measure the bulk 
delay for science targets. 

Figure [7] plots expected and measured ARV of GJ 699 over 11 nights in June, July, 
August and September of 2010, as well as the residuals. The radial velocities are expected 
to match the motion of the observatory with respect to the Solar System barycenter. The 



- 28 - 



fibers are treated independently, with independent template measurements. The template 
measurements, used as the zero point for radial velocity changes, were taken on June 17, 
2010 and have roughly the same signal-to-noise as the epoch measurements. Figure [H] plots 
the residuals for each fiber including the la formal uncertainties, and Figure |9] plots a 
histogram of residuals, normalized to the formal uncertainties. 

The root-mean-square (RMS) of the residuals is 36 m s~^ in Fiber A and 37 m s~^ in 
Fiber B. This is 1.76 times larger than the median formal uncertainty of 21 m s~^. After 
inflating the formal uncertainties of the measurements by 1.76, and calculating the one-sided 
Kolmogorov-Smirnov (KS) statistic between the residuals and a normal distribution of 
mean and standard deviation 1 (N[0,1]), we found the significance of the KS statistic 
to be 0.94. This indicates that the residuals are not significantly different than a normal 
distribution. The cumulative probability distributions of the normalized residuals and a 
N[0,1] distribution are plotted in Figure [TOl 

The individual ARV measurements can be averaged together to reduce the residuals. 
Figure im plots the RMS of the residuals versus averaging bin size. The reduction in RMS is 
not as fast as would be expected from white noise only. We model the reduction in RMS as 
a function of white noise aw and systematic noise ap, such that RMS = a/ crw^/Nbins + cf^- 
The systematic noise indicates the best precision achievable over these timescales which 
cannot be reduced by binning measurements. We calculate a best fit where aw = 33 m s~^, 
and (Ti? = 15 m s^^. 

The residuals in Figure M appear to have low frequency systematic fluctuations, 
characteristic of "red noise" or 1/f noise. The epoch measurements are unevenly sampled 
in time which makes it difficult to make an accurate power spectrum of the noise. 
Nevertheless, the noise can be analyzed versus observation number, rather than time. The 
Allan deviation, a a, measures the correlation between repeated measurements within a 
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data set (e.g. iThompson et al.ll200ll ). and is given by: 



" 2(iV + l-2A;) 2 " (15) 

where iV is the total number of measurements and Ui is the residual for measurement i. 
Figure WI\ plots the Allan deviation vs lag A;, in units of number of observation number, for 
the TEDI measurements of GJ 699. If the residuals were purely white, one would expect 
the Allan Deviation to decrease with the square root of the lag. The residuals follow a 
white noise contour until a lag with an Allan Deviation of 19 m s~^. 

The reduction in RMS with bin size and the Allan deviation calculation are both 
consistent with a noise floor of between 15 and 19 m s^^. It is advantageous to average 3-4 
measurements, but beyond that one is rapidly approaching the noise floor, and there are 
diminishing returns. In the next section we describe simulations of TEDI data and analysis, 
and find that incomplete telluric calibration is likely the source of the noise floor. 



6. Telluric Calibration Errors 

We expect the largest source of error in the measurements to be due to insufficient 
removal of telluric contamination in the complex visibilities. TEDI does not fully resolve the 
narrowest features of either the stellar or telluric lines. When unresolved, high resolution 
features of the stellar and telluric spectrum will mix and alias to lower resolution features 
which will contaminate the complex visibilities. The effect of the telluric and stellar mixing 
on the ARV measurement will depend on stellar type, projected rotational velocity, radial 
velocity, airmass of the target, and the telluric line strengths. All of these parameters will 
change the aliased spectrum and will contaminate the complex visibilities differently. 

To investigate this effect, we simulated mixed stellar and telluric spectra and telluric 
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calibration spectra, and carry them them through the analysis procedure. We used a 
high-resolution model spectrum of a 3200 K main sequence st ar provided by Tray is Ba rman, 



computed from the PHOENIX model atmosphere code (e.g. 



Fuhrmeister et al. 



20051 ) ■ and 



a high-resolution model of the telluric tran smission c alculated for Palomar Observatory by 



Henry Roe using his custom BEATS code (jRoe 



2OO2I) 



We convolved the stellar model with a 1 km s~^ rotational broadening kernel before 
shifting it by an input radial velocity, multiplying it by the telluric transmission model, 
multiplying it by the 10 interferometer transmission combs corresponding to the 10 phase 
steps in a phase set, and finally convolving it with a Gaussian profile to simulate the 
TripleSpec line-spread function. We matched the signal-to-noise of the simulated spectra 
to that of the GJ 699 data, and added Poisson noise to simulate photon noise. After this, 
we fed the spectra into the current data reduction algorithm described in Section HI using 
the same wavelength regions of 1.48 to 2.15 fim for measuring the change in radial velocity. 
The simulated data does not contain fluctuations due to telescope pointing and guiding, 
nor TliAr lines, and the delays are assumed to be known exactly. 

The expected radial velocities of the GJ 699 measurements were used as inputs, 
including the bulk radial velocity of t he star relative to the Solar System barycenter of 



-106.8 km s ^, previously measured by 



EvansI (1l996l ). Eigure fT3] plots the resulting residuals 



and la uncertainties based on photon noise for simulations with telluric contamination 
and calibration, and for simulations without any telluric contamination or calibration. The 
RMS of the simulated residuals with telluric effects is 43 m s~^, roughly a factor of 2 higher 
than that expected from photon noise, whereas the RMS of the simulated residuals without 
telluric effects is 13 m s~^, matching that expected from photon noise. The residuals show 
slow fluctuations, which correlate strongly with the expected radial velocities from the 
motion of the telescope relative to the Solar System barycenter. We interpret these errors 
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as due to uncalibrated mixing of stellar and telluric lines. It is not certain whether this is 
the source of the systematic noise in the TEDI measurements, since the stellar and telluric 
models are not perfect representations of GJ 699 nor the true telluric transmission, but this 
analysis suggests it is a plausible explanation. 

6.1. Fiber Illumination Errors 

TEDI does not currently have a mechanism for supplemental scrambling of the output 
of the fibers, as discussed in Section 13.31 This could result in slightly different cavity 
illumination between the ThAr calibration lamp and the starlight, and can potentially 
introduce errors into the measured ARV. The effect was tested by moving an artificial 
source on the input fiber tip. To do this, we built a fiber-fed telescope which we mounted 
on the top of TEDI to imitate a 1 arcsecond "star" at nearly the same etendue as the 200" 
telescope. We illuminated the fiber with a krypton gas discharge lamp, which has infrared 
emission lines in H and K bands. Moving the illumination of the TEDI input fiber by 0.5 
arcseconds was found to introduce a 45 m s^^ systematic offset in the measured ARV. This 
could potentially introduce stellar radial velocity offsets if the internal TEDI guider were 
guiding the starlight to the edge of the fiber core, rather than the center. However, when 
the fiber is illuminated 0.5 arcseconds off-axis, the images of the krypton emission lines on 
the spectrograph detector show a clear double-peaked slit profile, due to an annular shape of 
the output fiber illumination combined with slit-losses. The slit profiles of GJ 699 data do 
not show double-peaked slit profiles. Since the worst case scenario introduces error on the 
order of our current RMS performance, 45 m s~^ compared to 37 m s~^, and we do not see 
the double-peaked profile, we believe this effect is not currently limiting our performance. 
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7. Conclusions and Discussion 

We have demonstrated that the combination of a variable delay interferometer and a 
resolution 2700 near-infrared spectrograph can achieve better than 37 m s"^ of precision 
with 5 minutes of total integration time per observation on a nearby M dwarf, and that 
the largest source of error is hkely insufficient cahbration of telluric hues. At the current 
performance level, TEDI is not strongly competitive with current visible-wavelength 
Doppler surveys of M dwarfs. Nevertheless, the demonstration of better than 40 m s~^ 
of precision represents an important milestone in the effort to achieve precise radial 
velocities at infrared wavelengths, especially considering that the spectrograph involved has 
a resolution of only 2700. In order to be competitive with visible- wavelength radial velocity 
surveys, this technique will need to overcome the effects of insufficient telluric calibration 
and be implemented with significantly higher throughput. 

7.1. Reducing Telluric Effects 

The effects of insufficient telluric calibration can be addressed by several means. One 
solution is to forward model the complex visibilities using high resolution stellar and telluric 
spectra, either measured by high-resolution spectrographs or modeled by radiative transfer 
code. Unfortunately, accurate, high-resolution models of late type stellar spectra are 
challenging to produce because of the complex molecular opacities which vary significantly 
with temperature, pressure, and thus photospheric depth. While telluric transmission 
models have improved greatly over the last several years, it is unclear whether they are 
accurate enough for a forward-modeling approach to precise infrared radial velocimetry. 

Another method of reducing the effect of insufficient telluric calibration is to construct 
high-resolution spectra by combining many bulk delays from the interferometer. By 
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combining the complex visibilities from multiple delays, one can re construct the initia l high 



20031 ) ■ In 



resolution spectrum, similar to a Fourier transform spectrograph (lErskine et al.l 
principal, this method can achieve arbitrary resolution at the expense of survey efficiency, 
since each delay requires a phase set. In the case of TEDI, it would require 30 phase sets 
at different delays to achieve an equivalent resolution of 100000, which is the resolution 
of CRIRES. With the stellar and telluric lines fully resolved, or nearly fully resolved, the 
effects of mixing could be calibrated and reduced. Although, at this point it isn't clear what 
resolution is required. 



7.2. Increasing Throughput 

As stated in Section [372| the peak throughput of TEDI is 1.5% in H band. Most of the 
losses occur in the interferometer optical relay system and are due to the sheer number of 
optical elements, the use of 50% of the beam for throughput monitoring, and focal-ratio 
degradation in the fibers. The product of these losses results in 5% peak throughput 
through the interferometer alone. An optimum implementation of the interferometer could 
have dramatically higher throughput than TEDI. A fully fiber-coupled interferometer, 
placed away from the telescope, would significantly reduce the number of optical elements 
required and be able to relay both interferometer outputs to the spectrograph. With fewer 
optical elements, and better coupling of starlight into fibers, increasing the throughput of 
the interferometer up to 25%, and the total throughput of TEDI to 7.5%, is realistic. 



7.3. Technique Development 



Externally dispersed interferometry at near-infrared wavelengths is a propitious 
technique for measuring precise radial velocities of M dwarfs, but significant challenges 



-34- 



remain. It is important to determine whether the effects of telluric contamination present 
a fundamental limit to the performance, or whether significant performance gains can 
be made with better telluric calibration. If the technique were to efficiently achieve 10 
m s~^ of RMS precision on an M dwarf, that could dramatically affect the direction of 
future exoplanet instrumentation. The scientific returns from precise near-infrared radial 
velocimetry are rich enough to warrant pursing the last factor of 4 in velocity precision 
needed to achieve this goal, and attaining this higher precision is promising given the clear 
avenues remaining for improvement. 
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Fig. 1. — Schematic and beam path of the TEDI Interferometer, described in detail in 
Section 13.11 The yellow shaded region corresponds to the beam path. The pick-off mirror 
at (1) directs starlight to a guider system which injects the starlight and ThAr calibration 
light into fibers at (7). The fibers lead to the interferometer at (8), and eventually to a 
mirror which redirects the beam down to TripleSpec at (19). This entire assembly sits above 
TripleSpec at the Cassegrain focus and fits inside the 39" diameter hole in the 200" primary 
mirror. 
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Fig. 2. — An example TEDI exposure, flat-fielded and dark-subtracted, showing mixed 
starlight- ThAr lines on the upper fiber, and ThAr lines alone on the lower fiber, annotated to 
show wavelengths. Each fiber receives constant ThAr emission fight to measure the changes 
in delay introduced by stepping one of the cavity mirrors with a piezo actuator, with the 
starlight nodded between the two fibers. The starlight and ThAr light is cross-dispersed 
across 5 orders onto 2 quadrants of a Hawaii-2 HgCdTe detector. Background effects such 
as thermal emission and OH emission lines from the Earth's atmosphere are removed by 
subtraction after the spectra have been extracted. 
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Fig. 3. — A typical TEDI spectrum of mixed stellar (GJ 699) and TliAr light in units of 
signal-to-noise. Each of the 5 TripleSpec orders are plotted independently, and they overlap 
around 1.00 and 1.25 fim. Ten spectra, called a "phase set", are combined to form the 
complex visibility (see Figure H]). The extracted spectra are not corrected for instrument 
efficiency or telluric transmission, since the radial velocity signal is in the moire modulation 
of the spectra, not the overall shape. Telluric effects are calibrated out after modeling the 
moire fringes. The steep drop in signal-to-noise at 2.2 fim is due to the combination of higher 
background and lower transmissivity of the dichroic used for guiding (element 3 in Figure 
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Fig. 4. — TEDI spectra versus wavelength and At for a small region in H-band. Top: 
Mixed stellar (GJ 699) and ThAr light. The ThAr emission lines dominate the signal, but 
underneath them lie the stellar spectra. Middle: Separated ThAr lines used to accurately 
and precisely determine the values of tq and At. Bottom: The stellar spectra separated from 
the ThAr lines, after being normalized to remove effects from pointing/guiding fluctuations. 
The moire fringes are most evident between 1.57 and 1.58 /im and between 1.60 and 1.61 
//m, and these are due to two CO2 absorption bands in the Earth's atmosphere. The less 
distinct fringes between 1.585 and 1.60 are due to the stellar absorption lines, and the phase 
of the stellar moire fringes is highly sensitive to a Doppler shift. A^, Bj, and are fit to 
each wavelength channel to model the moire pattern. Regions with strong ThAr lines or a 
poor fit to the complex visibility have been blocked out. 
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Fig. 5. — Scatter plot of the phase differences for each wavelength channel between the ARV- 
shifted template complex visibility and an epoch complex visibility of GJ 699. For clarity, 
the phase differences have been converted to radial velocity differences using Equation |H1 
The scatter is roughly 1 km s~^, but the measured radial velocity depends on the weighted 
average of all of the wavelength channels. The clustering of points indicates spectral regions 
with narrow stellar features and strong radial velocity signal, with 1.7 to 1.8 fim containing 
the most signal in this bandpass. We currently limit the bandpass to 1.48 to 2.15 fim. 
The two gaps in coverage are from excessive telluric contamination near 1.62 fim, which is 
intentionally ignored, and a gap in TripleSpec's wavelength coverage around 1.95 ^m. The 
measured ARV is that which shifts the template to minimize the weighted value of these 
residuals. 
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Fig. 6. — Top: Histogram of the residuals in Figure [5] normalized to their la uncertainties, 
with a normal distribution of mean and standard deviation 1 included (N[0,1]). Bottom: 
Cumulative distribution function (CDF) of the residuals, including that expected from the 
normal distribution. The difference indicates significantly more outliers than expected from 
a Gaussian distribution. We remove those wavelength channels with the highest 3% of 
residuals and repeat the ARV fit to reduce the effect of outliers. 
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Fig. 7.— Top: Expected and measured ARV's of GJ 699 with TEDI. The expected ARV's 
are calculated as the motion of the telescope relative to the Solar System barycenter. This 
includes a 30 km s~^ semi-amplitude component from the Earth's orbit about the Sun, and 
a 300 m s~^ semi- amplitude component from the Earth's rotation. Top insert: Detail of 4 
nights in June, indicating clear recovery of the Earth's rotation. Bottom: Residuals including 
formal Icr uncertainties. 
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Fig. 8. — Residuals for each TEDI ARV measurement of GJ 699 with formal la uncertainties, 
equally spaced in order of observation, and grouped by observing run. The median la 
uncertainty of the measurements is 21 m s~^. Top: Residuals for Fiber A, with a RMS 
variation of 36 m s~^. Bot: Residuals for Fiber B, with a RMS variation of 37 m s~^. The 
fibers are treated independently, each with its own template measurement taken during the 
June run. 
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Fig. 9. — Distribution of the ARV residuals for GJ 699, normalized by their formal la- 
uncertainties. The median formal uncertainty of the measurements is 21 m s"^, however the 
RMS is 37 m s~^. Normal distributions of mean and standard deviations 1 and 1.76 are 
also shown. We include a Normal distribution with a standard deviation of 1.76 to account 
for the difference between the RMS and median formal uncertainty, and this distribution 
shows significantly better correspondence with the data. 
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Fig. 10. — Cumulative probability distribution (CPD) of the ARV residuals for GJ 699, 
normalized to the formal errors multiplied by 1.76, which is the ratio of the RMS of the 
residuals to the median formal uncertainty. This accounts for the difference between the 
median formal error and the measured RMS of 21 m s~^ and 37 m s~^, respectively. The 
CPD of a Normal distribution with mean and standard deviation of 1 (N[0,1]) is also 
shown. The one-sided Kolmogorov-Smirnov statistic is shown, and has a corresponding 
probability of 0.94, indicating that the distribution of residuals does not differ significantly 
from a Normal distribution. 
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Fig. 11. — RMS of the ARV residuals for GJ 699 as a function of number of measurements 
binned together. Included are plots indicating what is expected from a white noise component 
aw and a systematic noise component ap- RMS(A^Bins)=-\/'^w'/'^Bins + where A^Bins is 
the bin size. The solid line indicates that which would be expected from only white noise 
matched to the RMS when A^Bins = 1- The dashed line is fit to the RMS(A^Bins) and indicates 
a white noise component 33 m s~^ which can be reduced by binning, and a systematic noise 
component of 15 m which cannot. 
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Fig. 12. — Allan Deviation for TEDI residuals on G J 699 for each fiber. TEDI measurements 
are not equally spaced in time, so the lag used in the calculation is the epoch number, 
as plotted in Figure [81 Included is the expected contour for purely white noise with a 
standard deviation of 33 m s~^. The residuals are consistent with white noise until a lag 
of 4, corresponding to an Allan Deviation of 19 m s^^. This suggests that the non-white 
contributions to the noise amount to 19 m s~^ level, roughly consistent with ap from Figure 
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Fig. 13. — Simulated TEDI measurements of GJ 699 including calculated la uncertainties 
based on photon noise. Top: Simulated residuals including telluric contamination and cali- 
bration, versus epoch number, showing an RMS of 43 m s~^ despite a median Icr uncertainty 
of 25 m s~^. Bottom: The same, but without telluric contamination or calibration, and the 
RMS matches that expected from photon noise of 13 m s~^. The simulations suggest that 
insufficient calibration of the telluric contamination is currently limiting TEDI performance. 
The la uncertainties are significantly smaller for the simulations without telluric contami- 
nation because the telluric transmission reduces the photon counts at the edges of H and K 
bands. 
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